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A bisection algorithm for the numerical Mountain Pass 
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Abstract. We propose a constructive proof for the Ambrosetti-Rabinowitz 
Mountain Pass Theorem providing an algorithm, based on a bisection method, 
for its implementation. The efficiency of our algorithm, particularly suitable 
, for problems in high dimensions, consists in the low number of flow lines to 

^ ) ■ be computed for its convergence; for this reason it improves the one currently 

used and proposed by Y.S. Choi and P.J. McKenna in 

C3 



Introduction 

' This paper deals with constructive methods to seek critical points which are not 

minimizers for functionals defined on Hilbert spaces. The existence of such critical 
points may be detected from the topological features of the sublevels, under some 
compactness conditions. The purpose of this work is to give a constructive version 
of the Ambrosetti-Rabinowitz Mountain Pass Theorem (see pQ), which is one of the 
most useful abstract critical point theorems and that has found relevant applications 
in solving nonlinear boundary value problems. Our constructive proof relies on 
■ an algorithm based on a bisection method; this algorithm, elaborated in different 

versions depending on the topology of the sublevels (see Algorithms [21 and 0J , can 

S be implemented numerically. In spite of its simplicity, the algorithm proposed is 
new and in terms of numerical optimization, it improves the one currently used and 
. proposed by Y.S. Choi and P.J. McKenna in [3] (see also U). The efficiency of our 

' algorithm consists in the low number of steepest descent flow lines to be determined 

. to obtain a good approximation of a critical point of mountain pass type; this fact 

makes the method particularly fit for problems in high dimensions such as those 
coming from the discretization of nonlinear boundary value problems in infinite 
dimensional spaces. Moreover, compared with the usual Newton's method, our 
algorithm has two main advantages: it always converges and it does not need an a 
priori good initial guess. 

For some numerical applications of Algorithms |21 and 01 we refer the reader 
to Chapter 2 of [2], where such algorithms are applied to the action functional 
associated to the n-body problem with simple or double choreography constraint 
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respectively. The theory exposed in this paper, permits the determination of a new 
solution for the 3-body problem in a rotating frame with angular velocity u) = 1.5. 

1. An iterative algorithm for critical points 

1.1. The steepest descent flow. We consider a Hilbert space X and a func- 
tional / : X — > R of class C 2 . Fixed c e I, we define the c-sublevel of / as the 
open set 

f c {x S X : f{x) < c} 

and the set of critical points of / as 

Crit(f) := {x G X : V/(af) = 0}. 

The point xq G C'rit(f) is a ZocaZ minimizer for the functional / if there exists 
7' > 0, such that f(x) > f(xo), Vx G _B r (iro); xo is a strict local minimizer if there 
exists ro > such that for every r < ro, inf f(x) > f(xo). 

x£dB r (xo) 

Let r\ : R+ x X — ► X be the steepest descent flow associated with the functional 
/ defined as the solution of the Cauchy problem 

-n{t,x) - 



(1.1) { dt IK ^> l + ||V/(r,(i,x))|| 

77(0, x) = x 

We say that a subset X C X is positively invariant for the flow r\ if {w{t, Xq), t > 
0} C X , for every x n S X . We term u-limit of x G X for the flow 77, the closed 
positively invariant set 

lo x = { lim n(t n ,x) : (t n ) n C M+ 

We now state two useful preliminary Lemmata concerning some properties of the 
steepest descent flow defined in <|1.1[1 : for their simple proofs we refer, if necessary, 

to 

Lemma 1.1. Let n be the steepest descent flow defined in ]!.!)) : let x G X, 
7 G (0, 1] and T > 0, then 

\{t G [0, T] : \\Vmt,x))\\ > 7 }| < /(*Wfa(?>)) , 

Lemma 1.2. Let r\ be the steepest descent flow defined in M.l)) : let x G X, then 

(i) u> x G Crit(f); 

moreover if c x := lim t ^ +oc f(i](t, x)) > —00, then 

(ii) 3(t n ) n ,t n — > +00 such that 

lim f(r)(t n ,x)) = c x and lim V f(n(t n , x)) = 0. 

n — >+oc n — >+oo 

1.2. Disconnected sublevels. Let c G R be such that the sublevel f c is 
disconnected, we term (F?)i its disjoint connected components 

f c = \jF?, FfnF? = <n, yi^j. 

i 

For every index z, we consider the basin of attraction of the set F? 

Ti := \x G X : u> x C F t c } . 
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Proposition 1.3. Let f c be a disconnected sublevel for the functional f. Let 
Ff be the disjoint connected components of f c and Ff their basins of attraction. 
Then, for every index i, the following assertions hold: 

(i) Ff is an open set; 

(ii) <9Ff is a positively invariant set; 

(iii) mf f( x ) > c; 

(iv) minimizers of f in dTf are critical points for f ; 

(v) critical points for f in the set dTf are not strict local minimizers. 

Proof. We prove properties (i)-(v) when i = 1. 

(i) Let x £ Ff , then uo x C Ff and, since 77 is a gradient flow and Ff is a 
connected component of the sublevel / c , there exists T such that r](t,x) £ Ff, for 
every t>T. The continuity of the flow 77 ensures that there exists S > such that 
7](T,B 5 (x)) C Ff and, hence, B s (x) C Ff. 

(ii) By the sake of contradiction, suppose there exist x £ (9Ff and T > such 
that := f](T,x) £ <9Ff; if it € Ff, then there exists f such that n(T,x) £ Ff, 
hence a; G Ff , but this contradicts (i). If 17 £ -X"\Ff , then there exists e > such 
that B e (xx) C -X"\Ff; hence there exists (5 > such that rj(T, Bs(x)) C B € (xt) in 
contradiction with the definition of the set Ff . 

(iii) Follows from the definition of the set Ff . 

(iv) Follows from (ii) and ui x C C'rit(f), for every x £ dTf . 

(v) Let x £ 9Ff be a critical point for the functional /; by the sake of contra- 
diction, suppose that a; is a strict local minimizer, hence x ^ 9Ff. Let r > be 
such that B r (x) n Ff = and a r > such that 

/(a;) > a r + /(a;), Vx £ dB r {x). 

By the continuity of /, the set 

U := {x £ B r (x) : f(x) < a r + f{x)}, 

is a neighborhood of x, hence there exists a point xjj £ U D Ff ; since the flow 77 is 
dissipativc, and f{dB r (x)) > f(xu), we claim that w xu C B r (x), which contradicts 
the definition of Ff . □ 

We term path a continuous function 7 : [0, 1] — > X. Given a pair of points 
xi, a; 2 <E X, xi 7^ x 2 , we define the set of paths joining x\ to x 2 as 

(1.2) r xi , X2 :={ 7 eC([0,l],I):7(0) = x ll7 (l)=i2}. 

Theorem 1.4. Lef / c &e a disconnected sublevel for the functional f. Let Ff 
be the disjoint connected components of f c and Ff their basins of attraction. Let 
Xi £ Ff , i = 1,2, and 7 £ T Xl , X2 ; then there exists x £ 7QO, 1]) H <9Ff . 

Proof. The first step is the description of an algorithm that, given a path 7 
in the set r xitX2 , selects a point x £ 7QO, 1]) n <9Ff . 

Algorithm 1. 

StepO. s? = 0, *° = 1. ^ = £l±il, x o =!Bl> x o =X2> x O m=l{s oj 

Step i . if u x i-i C Ff° , sf = s 1 - 1 , s\ = 4- 1 
else s\ = s\ , s\ = s\ n 1 

gi _|_ gi 

x\ = 7(4), x\ = 7(4), s l m = 1 - 2 , a4 = 7(s^) 
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We have then defined two sequences (s\)i, (s|)j such that 

o = s? < s\ < . . . < s\ < 4 < • • • < 4 < 4 = i 

and 

l-i - 41 = ^ ^5° o. 

Since (s\)i, (s l 2 )i are bounded monotone sequences 
Let x = j(s), hence 



lim s\ — lim s\ — s £ [0, 1]. 

-►4-oo i — >+oc 



lim x\ = lim x\ — x, 

i — > + oo i — > + oo 

and, necessarily x lies on the positively invariant set dTy . □ 

Corollary 1.5. In the same conditions of Theorem \l.J\ let x € 7([0, l])n9J-"f , 
f/ien > c and £/iere exists a sequence (x n ) n = r](t n ,x) C ftFf smc/i that 

lim V/(x„) = 0, lim f(x n ) = 

n — >+oc n — > + oo 

Proof. Since dTy is a positively invariant closed set, w s C dTy and fiy>x) > c. 
We use Lemma fl .21 to define the sequence {x n ) n . □ 

Unfortunately, the proof of Corollarv ll.5l is not constructive in the sense that 
it does not provide a method to determine the critical set uj s . The reasons why 
we can not have an implementable proof of this result are, first, that we can not 
determine precisely, in a finite number of steps, the point a;, since it is the limit of 
the sequence (x 1 l ) n in Algorithm ^ Second, w 2 is defined as a limit for t — > +oo 
and we are not able to determine the value /(wj) numerically. 

Although the following result is a consequence of Corollary 11.51 its relevance 
consists in its constructive proof, that gives a method to determine a critical point 
for the functional / at a level higher than c. 

Corollary 1.6. In the same conditions of Theorem \l.J\ let x £ 7([0, V\)C\dT1, 
then there exists a sequence {y n )n C X, y n :— r](T n ,Xi), such that 

lim Vf{y n ) = 0, and c < f(y n ) < f(x), Vn € N. 

n — >+oc 

Proof. Let (x™)„ be the sequence defined in Algorithm ^ { x i)n C T{ n 
7([0, 1]) and converging to x £ dT{. Since [x\) n C Ty, we can define 

T n :=inf{t>0:/(»/(i,a;?)) <c} 

and 

7n := inf ||V/(77(*,a?))||. 

te[o,T„] 

Hence, using Lemma 1 1.1 1 (we can suppose 7„ < 1) we deduce 

T n = \{t £ [0,T B ] : \\Vf( V (t,x^)\\ > 7 „}| < - C) 

In 

and we can conclude that 

l2(f(x?)-c) 



In < 
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Since (x") n — ► x as n — > +00 and f(u>s) > c, we have 

lim T„ = +00 and lim 7„ = 0. 

n — > + oo n — >+oo 

We deduce the existence of a sequence (T n ) n such that (T n )„ £ [0, T n ] and, defining 
2/n := v(T n ,Xi), we have 

lim ||V/(y„)|| = 0. 

n— »+oc 

□ 

The proof of Corollary 11.61 shows that Algorithm ^ can be improved in the 
following Algorithm El to obtain the sequence (y n )n- 

Algorithm 2. 

Step 0. s° = 0, s° = 1, .s°„ = ^±£2 , =x u x° 2 =x 2 , x° m = 7 («°J 
Step i . if t^.-i C , s\ = s^ 1 , s| = 4 _1 
else s\ = , S2 = s m 1 

5* -)- 

xl=7(4). 4=7(4). 4" 1 2 2 . 4»=7(4i) 
1} := inf £>():/(»#, a:!)) <c} 
Tj := t e [0,Ti] that minimizes ||V/(?j(t, x{))\\ 
Vi := T](fi,x\). 

The sequence (y n )n defined in Algorithm [3 will converge, once we impose some 
additional compactness conditions on the functional /. In this sense we give the 
following definitions. 

Definition 1.7. A sequence (x m ) m C X is termed a Palais-Smale sequence in 
the interval [a, b] for the functional f if 

«</M<6,VmeN and Vf(x m ) m ^X°° 0. 

The functional / satisfies the Palais-Smale condition in the interval [a, b] if ev- 
ery Palais-Smale sequence in the interval [a, b] for the functional /, (x m ) m , has a 
converging subsequence x mk -> 10 G X. Similarly, a sequence {x m ) m C I is a 
Palais-Smale sequence at level c for the functional f if 

j(a: m ) — ► c and \/J(x m ) > 0. 

The functional / satisfies the Palais-Smale condition at level c, (PS) C , if every 
Palais-Smale sequence at level c for / has a converging subsequence. 

Remark 1.8. Corollary 11.51 ensures the existence of a Palais-Smale sequence 
at level /(w 2 ) for the functional /. When the functional / satisfies the (PSW^-j, 
we conclude that there exists a critical point x for / such that f(x) = /(wj). 
Corollary 11.61 implies the existence of a Palais-Smale sequence for / in the interval 
[c, f(x)]. When the functional verifies the Palais-Smale condition in [c,f(x)], it 
ensures the convergence of the sequence (y n ) n constructed in Algorithm^ 

Figuren]shows the implementation of Algorithmic In the first graph, each line 
represents the values of the action functional on the steepest descent flow starting 
at the first elements of the sequence (x\)i. Remark that, as the index i increases, 
the interval of time in which the action function is approximately 5.2 increases, too. 
This is due to the fact that the line of the flow departing from x\ passes each step 
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closer to the desired critical point, whose action level is approximately 5.2. In the 
second picture, we see the value of the norm of the gradient of the action functional 
on the same curves. This value decreases till the point yi, then it increases when 
the line r](t, x\) departs from a neighborhood of the mountain pass point, to reach 
a local minimizer. Remark that both figures represents just the first 1000 steps in 
the t- variable. 

We now point out some cautions to be taken in the implementation of Algorithm 
13 After a certain number of steps, depending on the distance between x\ and x%, 
the points x\ and x\ may not be any more numerically distinct; moreover, taking 
into account the numerical errors in the integration method and a possible strong 
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sharpness of the graph of /, the lines of the steepest descent flow starting from x\ 
and x\ can become undistinguished. Hence, each time wc implement Algorithm [3 
we a priori fix a maximum number of iterations, N max , to avoid such numerical 
obstructions and to make the distance dist(x^ Vmax , x^™"*) <C e, where e is a fixed 
small positive number. Fixed e,7 > we propose the following algorithm that 
allows us to approach a locally optimal path joining the starting points xi and X2. 

Algorithm 3. 
Step N^ ax — N max ; 

Algorithm 2(xi,X2, N^ ax ; Vn^J 
Step k if \\^f(y N ^)\\ <7> STOP 

f / TV'" -1 iV* _1 \ 1 

else T def := inf U > : dist \ r](t 7 x 1 ma *),r\{t 1 x 2 max ) j > ^fer j '• 

Xj :=-q(T def ,x^ max ), J = 1,2; 
N k ■= N k ~ x - 1 • 

max ' max > 

Algorithm 2(xi, x 2 , N^ ax ; Vn^J 

It must be said that the choice of the second connected component where Xi 
lies is forced by the problem itself. Nevertheless, we can connect X\ to any given 
connected component by juxtaposition of a finite number of locally optimal paths. 

1.3. Sublevels with non trivial fundamental group. We recall that a 
path 7 is a continuous function 7 : [0, 1] — > A. Given a path 7, we define the path 
—7 as 

-7(«) := 7(1 - a), Vs G [0, 1] 
and for every pair of paths 71,72 such that 7i(l) = 72(0) we define the path 
7 = 7j o 72 juxtaposition of 71 , 72 as 

7{S) (7io l2 )(s) . I ^-i) SG [1,1] 

A path 7 : [0, 1] — > X is a loop if 7(0) = 7(1). We say that a loop 7 is contractible 
in X if there exist x £ X and a continuous function h : [0, 1] x [0, 1] — > X such 
that /i(0,s) = 7(a), /i(l,s) = x, Vs e [0,1], and h(X,0) = h(X,l), VA G (0,1). We 
observe that, for every A G [0, 1], h\(-) :— h(X, •) is a loop in X. 

Let Y be a subset of A, we say that Y is simply connected if every loop in Y 
is contractible. 

We consider a functional / : X — ► ffi, satisfying the following hypotheses 
(hi) / is bounded below, suppose f(x) > 0, Vi e A; 
(h2) there exists c such that the sublevel f c ^ is not simply connected; 
(h3) / verifies the (PS) Co for every < c < c. 

In the sequel, we will work in a connected component of the sublevel / c ; we 
will still term it f c . 

Let x € f c be a strict local minimizer for / and r > such that the ball B r (x) is 
contractible in f c , i.e. for every x G B r (x) the set {(1 — A)x + Ax : A G [0, 1]} C f c . 
We consider the sublevel / c+£ , where c = /(x) and e > 0, such that -F| +e C B r (x), 
where F£ +e is the connected component of f c+e containing x. We observe that 
i^ +£ is contractible in f c . 

For every x± G / c such that w I1 = x we define the instant 

T X1 := mi{t > : r){t, x x ) G F| +£ } > 0. 
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We define the following paths 

(1.3) 7^! (A) := rj(T Xl X, xi), J 2 X1 (X) := (1 - X)r,(T Xl , Xl ) + Xx, 
Definition 1.9. We term descending path associated to x\ the path 

<*X1 ■= (7x! °7xJ 

where j xi ,7™ are defined in (|1.3(l . 

Remark 1.10. We observe that, defining X := {x £ X : uj x = x} and T x x := 
{7 G L^s : x e X}, the maps 

X —>■ R_|_, a;i 1— > T Xl and X — > T x x , x i l— * 
are continuous on the set X. 

Definition 1.11. Let 7 be a path in f c such that w 7 (o) = w 7 (i) — 2;; we define 
the loop 7^, the descending loop associated to 7, as 

(1.4) 7„ := [(-a T (o))°7]°a T (i). 

Definition 1.12. We say that a path 7 in f c is rj- contractible in f° if 

(a) w 7 ( ) = w 7 (i) = x; 

(b) the loop 7,, is contractible in / c . 

Lemma 1.13. Let x £ f c be a strict local minimizer for the functional f and 7 
a path in f € . If to^u) = x, Vs £ [0, 1], then 7 is 77- contractible in f c . 

Proof. To prove that the path 7^ defined in (|1.4fl is contractible, we define a 
continuous function h : [0, 1] x [0, 1] — > f c such that /i(0, •) = 7 r) (-) and h(l, ■) = x. 
For every A € [5,1] we define the path 

/3 A ( S ):=a 7(1) (2(l-A) S ), s e [0, 1] , 

hence, the following function h satisfies our requirements 



□ 

Theorem 1.14. Let f be a C 2 functional on a Hilbert space X. Suppose that 
f satisfies conditions (hi), (h2), (h3). Then one of the following situations occurs: 

(i) f c has a continuum of global minimizers; 

(ii) there exists x € Crit(f) n f c that is not a strict minimizer. 

Proof. When the set Crit(f) n f c is not a continuum of global minimizers, 
hypotheses (hi) and (h3) imply the existence of at least one strict minimizer ief 
for the functional /. Let 7 C f c be a loop such that 7(0) = 7(1) = x and suppose 
that 7 is not contractible in f c . Hence, since 7, ; (s) = 7(5) for every s £ [0, 1], 7 
is not ry-contractible in f c . Using Lemma f 1.1 31 we can conclude that there exists 
Xy G 7([0, 1]) such that lu x 7^ x. The critical point uj x could be a strict minimizer 
or not. In the first case, we use Algorithm |21 to prove the existence of a critical 
point for / that is not a strict local minimizer; otherwise oj x ^ is the searched critical 
point for /. 

We now show an algorithm that determine a point x-y £ 7([0, 1]) such that 
uj Xi 7^ x. Let 7 C f c be a not contractible loop such that 7(0) = 7(1) = x. We 
have already remarked that 7 is ?/-contractible. 



f [(-a 7 (2A)) °7]° a 7 (i)j A £ 
\ (-/3a)°/3a, Ae 
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Algorithm 4. 

Step 0. 7o := 7. = 7o (|) 
if uj Xq ^x, STOP" 

if 70 ([0,5]) is 77-contractible, 71 := 70 ([§, l]) 
else 71 := 70 ([0, |]) 

Step i . Xj = 7i (i) 

if ^5, STOP 

if 7i([0)5]) is ?7-contractible, 7 i+ i := 7, l]) 
else 7i + i := 7^ ([0, 5]) . 

Using Algorithm 01 we can directly find a point whose w-limit is not x, in this 
case the i-loop is stopped. Otherwise, we define a sequence of paths that are 

not ^-contractible. In this case, let x® — 7i(0) and x\ — 7i(0) be the initial and 
end points of these paths, then the sequences (x®)i, {x\)i C 7([0, 1]) and 

lira x\ = lim x\ — S. 

i — > + oo i — > + oo 

Arguing as in the proof of the convergence of Algorithm^ we deduce that necessar- 
ily Log x; if cjg is a strict minimizer, we use Algorithm El to deduce the existence 
of a critical point that is not a strict local minimizer. □ 



2. The Mountain Pass Theorem 

Our goal now is to prove a version of the Mountain Pass Theorem (for a detailed 
theory on this subject we refer to [BE]). ^ Mountain Pass Theorem concerns itself 
with proving the existence of critical points which are not strict local minimizers 
for the functional /; using Corollary II .51 and Proposition 11.31 we are now able to 
prove the following 

Theorem 2.1 (Mountain Pass Theorem). Let f be aC 2 functional on a Hilbert 
space X. Let X\,X2 S X, let T Xl >X2 be the set of paths defined in H1.2(l and Co the 
level 

(2.1) c := inf sup /( 7 ( a )). 

76r xi , X2 se [ 0j l] 

such that 

(2.2) co >max{/(zi),/(z 2 )}. 

// the functional f satisfies the Palais-Smale condition at level Co then there exists 
a critical point for the functional f at level Co , that is not a local minimizer. 

PROOF. The definition l|2.1(l of the level Co and condition 12.2J1 . imply, first, 
that the sublevel f c ° is disconnected, second that for every k 6 N, there exist 
7fc € T XuX2 such that 

sup/(7fe(s)) < c + -. 

From Corollary 11.51 we deduce the existence of an element Xk € 7/c([0, 1]) H dT'y 
such that 

co < f(ux k ) < c + 7. 
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Let (x" k )n C T{° converging to Xk] following the proof of Corollary 1 1.61 we deduce 
the existence of an index mef) and an instant T Uk such that, y Uk :— rj(T nk , x" k ), 
where x" fc € { x ik)n, verifies 

|V/(y„J|| < ^ and c < /(y„J < c + -. 

Since / satisfies the Palais- Smale condition at level cq, we deduce that there exists 
a critical point, y, for / at level Co, limit of the sequence {y nk )k as k tends to +00. 

We conclude the proof showing that there exists a sequence (y„ k )k, converging 
to a critical point y* that is not a minimizer and such that /(j/^ fe ) — * f(y*) = Co, 
with f(Vn k ) = c o - \ < c - Let 

T* :=iDf(t>0:/fa(i, *?*))<<*- 1 



k 

hence, defining y* lk := r q(T* k ,x 7 l h ), we have f(y nk ) = Co — ^. If we can prove that 

lim v/(y;j = o, 

using the Palais-Smale condition at level c for /, we conclude the proof. By the 
sake of contradiction suppose that there exists e > such that || V/(y* )|| > eo, for 
every k € N. Hence, since the Palais-Smale condition at level Co holds, there exists 
S > such that \T* - T„J > 5 and \\Vf(r)(T* k - t,< fc ))|| > e' > 0, Vi <E [0,6). 
Then, using Lemma ITTTI 



-ti" ,:ci ))l1 - v 1 2 



/ 2(/(y nt )-co + l) < 1 



<$ 2 " 

and the right-hand side tends to as A; tends to +00. □ 
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